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We study heat conduction in (n, 0) / (2n, 0) intramolecular junctions by using molecular dynamics 
method. It is found that the heat conduction is asymmetric, namely, heat transports preferably 
in one direction. This phenomenon is also called thermal rectification. The rectification is weakly 
dependent on the detailed structure of connection part, but is strongly dependent on the temper- 
ature gradient. We also study the effect of the tube radius and intramolecular junction length on 
the rectification. Our study shows that the tensile stress can increase rectification. The physical 
mechanism of the rectification is explained. 

PACS numbers: 66.70.+f, 44.10.-K, 61.46.Fg, 65.80.+n 



I. INTRODUCTION 

In past two decades, the study of heat conduction in 
low dimensional model systems has enriched our under- 
standing on heat conduction from microscopic point of 
view 1 -. In turn, the study has also lead to some inter- 
esting inventions for heat control and management de- 
vices. For example, the heat conduction in nonlinear lat- 
tice models 2 -' 3 -' 4 - demonstrates rectification phenomenon, 
namely, heat flux can flow preferably in one direction. 
Furthermore, the negative differential thermal resistance 
is also found and based on which, a thermal transis- 
tor model has been constructed* 6 - Most recently, the 
two segment model of thermal rectifier proposed in Ref 
EH has been experimentally realized by using gradual 
mass-loaded carbon and boron nitride nanotubespl These 
works have opened a new era for heat management and 
heat control in microscopic level. 

On the other hand, there have been increas- 
ing studies on heat conduction in real nano scale 
systems^. For example, the thermal conduction 
of single walled carbon nanotubes (SWCNTs) has 
attracted both theoretical 9 ' 1 ^ 1 ' 12 ' 13 ' 14 ! 15 ' 16 ! 17 ' 18 and 
experimenta l 19 i 20 ' 21 i 22 ' 23 i 24 i 25 attentions. Almost all ex- 
periments and numerical simulations have payed their 
attention to the extremely high thermal conductivity of 
SWCNTs. The dependence of the thermal conductiv- 
ity on the length or chirality of SWCNTs is also exten- 
sively studied theoretically. Therefore, one may asks, 
how can we make heat controlling devices, i.e., the ther- 
mal diode or the thermal transistor, from SWCNTs and 
their derivatives? 

In the practical applications, SWCNTs may have 
many kinds of impurities such as element impuri- 
ties, isotopic impurities 26 , topological impurities, etc. 
Among these numerous derivatives of SWCNTs, the 



SWCNT intramolecular junctions (IMJs), which are 
formed by introducing the pentagon-heptagon rings in 
SWCNTs, have been expected to be an important 
structure in the applications. It is well-known that 
the electronic properties of the SWCNT IMJs have a 
close relationship with their geometrical or topological 
characteristicSi 27 ! 28 ' 2 ^' 3 ^' 3 ^' 3 ^' 3 ^ 3 ^ It is reasonable that 
the thermal transport behavior of SWCNT IMJs can 
also be altered by their geometrical characteristics. Com- 
pared with the electronic properties, the thermal trans- 
port behavior may have not so close relationship with 
the detailed local geometrical arrangements because it is 
mainly influenced by the long wavelength phonons. But 
the non-equilibrium thermal transport behavior of IMJs 
can be affected by the high-frequency optical phonon 
modest which can reflect the information of the lo- 
cal defects. So it is interesting to investigate the non- 
equilibrium thermal transport of different SWCNT IMJs. 

In this work, the thermal rectification in (n, 0)/(2n, 0) 
IMJs and its dependencies on tube radius, IMJ length, 
and external stress, e.g., tensile and torsional stress are 
studied. A complex structure of 'peapod' structure and 
IMJ is also studied to investigate the effect of periodical 
potential on the thermal rectification. 



The paper is organized as the follows. In Sec. II we 
introduce the basic structure of (n, 0)/(2n, 0) IMJ and 
our numerical method. In Sec. Ill, the main numerical 
results are discussed. In Sec IV, we show that the ex- 
ternal stress can improve the rectification. In the same 
section we also show the results of (n, 0)/Cgo@(2n, 0) 
structures, which can be regarded as a combination of 
peapod structures and (n, 0)/(2n, 0) IMJs. In Sec. V, 
we give some concluding remarks. 



2 



II. MODEL AND METHODOLOGY 

A typical (n, 0)/(2n, 0) IMJ structure is depicted in 
Fig. 1, in which the index n equals to 8. The struc- 
ture contains two parts, namely, a segment of (n, 0) 
SWCNT and a segment of (2n, 0) SWCNT. For sim- 
plicity, the lengths of the two segments are almost equal 
in our calculations. The two segments are connected by 
m pairs of pentagon-heptagon defects. Because (n, 0) 
and (2n, 0) tubes have the common rotational symme- 
try C271, thus when n can be divided exactly by m, we 
can adjust the defects to make the two segments connect 
straightly. Especially, if m = n, the pentagon-heptagon 
defects can be arranged in the connection part so that 
the rotational symmetry of the IMJ is C n . On the other 
hand, it is necessary to define a basic system length be- 
cause the thermal conductivity depends on the system 
length. 26 After setting the lengths of (n, 0) and (2n, 0) 
tubes to be the same, we can define a basic total length 
Lq when both of the two segments contain 24 periods. All 
the structures are fully optimized before further molecu- 
lar dynamic (MD) calculations, and then we can obtain 
Lq ps 20 nm. 
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FIG. 1: (Color online). A typical (n, 0)/(2n, 0) structure. 
Here, n = 8, and the number of pentagon-heptagon defects is 
m = 4. The regions marked as T are fixed in MD process. 
The regions marked as 'IF are put in the heat baths. 

In this kind of structure, the outmost one period of 
each heads (colored by blue and marked as region 'I' in 
Fig. 1) are fixed in MD process. Then two periods of 
each ends (illustrated by orange color and marked as re- 
gion 'II' in Fig. 1) are put in the heat baths, which are 
realized by the Nose-Hoover thermostat^ 6 - The temper- 
atures of the thermostats at left and right heads are Tl 
and Tr, respectively. For convenient, here we introduce 
two quantities: AT = Ti -~ Tr , and (T) = Tl 1 2 Tr . In this 
work, (T) is always kept at 290 K. 

The C-C bonding interactions are described by the 
second-generation reactive empirical bond order (REBO) 
potential^! which is the most recent version of the 
Tersoff-Brenner type potential, combining advantages of 
the two sets of parameters in the earlier version^ 

The velocity Verlet method is employed to integrate 
the equations of motion with the time step of 0.51 fs. 
The typical total MD process is 5xl0 6 steps, which is 
about 2.55 ns, and the statistic averages of interesting 
quantities start from half of the MD process, i.e., 2.5x 10 6 
steps are used to relax the system to a stationary state. 

The instant temperature of atoms i is defined as 

T * (*) = §t { Vx W 2 + y y (tf + v * W 2 ) ' where V W is the 

time-dependent velocity, m$ the mass and fcf, the Boltz- 
mann constant. This is a result of energy equipartition 
theorem. 



The thermal flux is obtained from the thermostats, i.e., 
the total work from the thermostats can be regarded as 
the heat flux runs from thermostats to the system. If the 
work is positive, then the heat flux is also positive. In the 
scheme of the Nose-Hoover thermostat r 6 - the equation of 
motion for the particle i in heat bath is: 
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where pi is the momentum and /j is the force applied on 
the atom, g is the number of degrees of freedom of the 
atoms in the thermostat. Q — gkbTr 2 , where r is the 
relaxation time. In our simulation, r is kept as 4 ps. The 
heat bath acts on the particle with a force — thus the 
power of heat bath is — £ , which can also be regarded 
as the heat flux from heat bath, i.e., 



Ji 



Pi - Pi 



rn, 



(2) 



The total heat flux from the heat bath to the system 
can be obtained by J = J^i^i) where the subscript i 
runs over all the particles in the thermostat. The final 
thermal conduction is k ■ s = 2 at/l = 5SF> where j is 
the heat flux density, s the area of cross-section, and L 
the system length. In fact this definition is equivalent 
to the usual definition: J = Jj (f) (t), (t) is the 

i 

instant total energy of particle i, but Eq. ([2]) is much 
simpler for computational simulation. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 



First of all, we investigate the effect of defects in con- 
nection region. The defects in connection region deter- 
mine the property of the interface between two tubes, 
and as a result, they will affect the heat transport prop- 
erty of the system. However, the question whether they 
can affect the rectification is still open. 

We study the (8, 0)/(16, 0) structures as examples. 

Generally speaking, they are many methods to connect 
an IMJ, thus we only consider two conditions here. The 
first one is the coaxal straight IMJs. The possible number 
m of pentagon-heptagon defect pairs for a coaxal IMJ is 
8, 4 or 2, but m = 2 makes the cross-section of the tube 
part highly deformed from circle, so we only consider m = 
4 and 8. The second condition is the simplest IMJ, i.e., 
two tubes are connected with each other by using only 
one pair of heptagon-pentagon defects (m = 1)^ The 
top views and side view of the three structures (m = 8, 
4 and 1) are plotted in Fig. 2. It can be found even in 
the simplest zigzag/zigzag IMJ (m = 1), there is still a 
small angle between the two segments. 

The total heat fluxes J versus the normalized temper- 
ature difference AT/(T) is plotted in Fig. 3a. Positive 
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FIG. 2: (Color online). The top views and side view of 
two (8, 0)/(16, 0) IMJs with different connection methods, 
(a) With 8 pairs of pentagon-heptagon defects, (b) With 4 
pairs of pentagon-heptagon defects, (c) With only one pair of 
pent agon- heptagon defects. 

J means the heat flux flows from (n, 0) tube to (2n, 0) 
tube. We do not show the heat flux density here, because 
it is difficult to define the area of the cross-section of the 
whole IMJ. 
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FIG. 3: (Color online). The thermal transports of (8, 
0)/(16, 0) SWCNT IMJs with different numbers of pentagon- 
heptagon defects in the connection region, (a) The heat flux 
versus temperature difference, (b) The thermal rectification 
versus temperature difference, m is the number of pentagon- 
heptagon defects. Error bars are also plotted. 

It can be found from Fig. 3a that the total heat flux 



for m = 4 is larger than that for to — 8. It is straight- 
forward because the pentagon-heptagon defects are topo- 
logical defects, they scatter the phonons and cause larger 
interfacial thermal resistance. As a result, more topolog- 
ical defects will reduce the thermal conductance and the 
heat flux. But the total heat flux for to = 1 is close to 
that for to = 4. This is due to the big difference of the 
connecting region of these IMJs. The length of the con- 
necting region in IMJ with m — 1 is much longer than 
those in the IMJs with m — 4 or 8. The long connecting 
region can cause big thermal resistance, as a result, the 
thermal resistance in to = 1 IMJ is not much smaller 
than that in to = 4 IMJ. 

When the temperature difference is very small, the 
heat flux only shows weak asymmetry. The asymme- 
try becomes larger when temperature difference is larger. 
This characteristic is seen clearer in Fig. 3b, in which the 
thermal rectification \J n /J P \ = —Jn/Jp is drawn. 

In Fig. 3b, J n means that the heat flux is of nega- 
tive sign, vice versa, J p means that the heat flux is of 
positive sign. It can be seen that, an increase in tem- 
perature gradient |AT| leads to an increase of thermal 
rectification. This is because the rectification is a result 
of non-equilibrium transport, which is mainly determined 
by optical phonons. Only when the temperature differ- 
ence is large enough, the optical phonons can be excited 
and contribute to the heat conduction. More interest- 
ingly, when the temperature difference is large enough 
(for | AT | w 0.8(T)), the rectification begins to decrease. 
But, the temperature difference is so large that this re- 
sult may not be hold in more restrict calculations due 
to lose of local thermal equilibrium and the quantum ef- 
fect in the lower temperature head. Another meaningful 
feature is that, although the absolute values of the heat 
fluxes are completely different, rectifications are almost 
the same in IMJs with m — 4 and 8, whose connecting 
region is short. Namely, the rectification is weakly de- 
pendent on the detailed structure of the interface when 
the connecting region is short enough. And the rectifi- 
cation in to = 1 IMJ is largest among the three kinds 
IMJs. 

In order to understand above rectification phe- 
nomenon, we consider the power spectra of the atoms 
around the connecting part. The atoms to be studied are 
illustrated in the Fig. 4a. They are selected to be next 
to the defects along the axes. The corresponding phonon 
spectra are presented in the Fig. 4b. And the overlap- 
ping of the phonon spectra of the two atoms around the 
connecting parts have been emphasized by shadows. 

It is found that when the direction of the tempera- 
ture gradient is exchanged, the area of the overlap region 
is also changed. More importantly, the overlap area is 
larger when AT/ (T) is negative. This result corresponds 
with the trend of the heat flux. So we can state that 
in this real system, the relationship between the overlap 
area and the absolute value of the heat flux is same as 
that in the one dimensional nonlinear lattice systems^i^, 
in which it is found that matching/mismatching of the en- 
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FIG. 4: (Color online), (a) The side view of the (8, 0)/(16, 
0) IMJ (m = 8), where the atoms whose power spectra are 
recorded are emphasized by balls, (b) The corresponding 
phonon power spectra. The overlapping of the phonon spec- 
tra of the two particles around the connecting parts have been 
emphasized by shadows and bordered by dark lines. 



ergy spectra around the interface is the underlying mech- 
anism of the rectification. 

On the other hand, the most obvious change appears at 
400-1000 cm" 1 and 1600-1800 cm" 1 . The m = 4 IMJ 
is also studied, and similar conclusion can be obtained. 
So the optical phonon modes are important for the rec- 
tification. In addition, in Ref. [35l , the authors show that 
in carbon nanotubes with finite length where the long- 
wavelength acoustic phonons behave ballistically, even 
optical phonons can play a major role in the non-Fourier 
heat conduction. The dispersion relations of the SWNT 
show that, in the intermediate range of the normalized 
wave vector 0.1 < k* < 0.9, some of the phonon branches, 
especially the ones with relatively low frequency, have 
group velocity comparable to the acoustic branches. 

Furthermore, the thermal conductivity can be ex- 
pressed as k — Y^, lxqC q (uj\) v\ q , where A is a set of quan- 

turn numbers specifying a phonon state, I is the phonon 
mean free path, v\ q is the magnitude of the phonon 
group velocity along the direction of the heat flow, and 

C q M = k B is the thermal ca- 

pacity of lattice wave with wave vector q and angular 
frequency u>. Thus, the contribution from the high fre- 



quency modes is weak because their group velocities are 
much smaller than those of low frequency modes. So the 
modes at 400—1000 cm -1 are more important than the 
modes at 1600-1800 cm" 1 . 

The frequencies of the low frequency modes have been 
found to be inversely proportional to the tube radius4^ 
Thus, the low frequency modes, e.g., the radial breath- 
ing mode, have different exciting temperatures in (16, 
0) and (8, 0) tubes. On the other words, their partici- 
pation in transport process occurs at different tempera- 
tures. This fact further explains why these low frequency 
optical phonon modes are important for the rectification. 

Based on above discussions, now we can understand 
the following phenomena. Firstly, the rectification almost 
does not change in the m — 4 and 8 IMJs. Secondly, the 
rectification in m = 1 IMJ is largest in the three IMJs. 

As has been mentioned, the modes at 400—1000 cm -1 
are most important for the thermal rectification. These 
modes have larger wave lengths than the high frequency 
modes. The connecting parts of the two IMJs with m = 4 
and 8 are short, as a result, the defects have weak effect 
on the effective optical modes. In Ref. |4lj, it is found that 
the high frequency vibrational modes caused by topo- 
logical defects are localized states. And the higher fre- 
quency localized modes show smaller spatial dimension, 
which means they can hardly entangle with the trans- 
port phonon modes. Thus in the short connecting IMJs 
(rn = 4 and 8), the arrangement of the detailed defects 
has weak influence on the rectification. In contrast, the 
connecting region is much longer in the m = 1 IMJ, and 
the distance between the pentagon and heptagon rings is 
relatively long, which means the connecting region can 
affect those middle frequency phonon modes, which have 
longer wave lengths. Therefore, the rectification is largest 
in the case of m = 1 IMJ. 

However, the rectification in current structures is still 
small, which might limit the application of the IMJ as a 
thermal rectifier. So we will try to find the factors which 
can help improve the rectification. 

Above all, we investigate the thermal conductance of 
the IMJs with the different index n numerically. 

In order to make the results comparable with each 
other, the number m of pentagon-heptagon defect pairs 
are set to be n. Considering the computational consump- 
tion, only n — 7, 8 and 9 are considered in Fig. 5. 

According to Fig. 5a, the total heat flux increases when 
n increases, but this does not mean that the thermal 
conductivity for big n is also large because it should be 
scaled by the area of cross-section. In fact, it has been 
shown that the thinner SWCNTs have higher thermal 
conductivity in literatures! 15 i 18 i 23 

Another interesting feature is that, when index n in- 
creases, the asymmetric heat flux ratio does not always 
increase. In Fig. 5b, the largest rectification appears 
when n = 8. So this reminds us that the heat flux rec- 
tification is not only induced by the radius difference of 
two segments. When tube radius increases to a suffi- 
cient large value, i.e., n is large, the vibrational density 



5 




0.4 0.6 

\AT\/{T) 

FIG. 5: (Color online). The thermal transports of (n, 0)/(2n, 
0) SWCNT IMJs, n = 7, 8 and 9. (a) The heat flux versus 
temperature difference, (b) The thermal rectification versus 
temperature difference. In this figure, m = n. Error bars are 
also plotted. 
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FIG. 6: (Color online). The thermal transports of (8, 0)/(16, 
0) SWCNT IMJs with different length, (a) The heat flux 
versus temperature difference, (b) The thermal rectification 
versus temperature difference, m — 4. Lq is defined in text. 
Error bars are also plotted. 



of states (VDOS) of carbon nanotube approaches to that 
of the graphite. But when n is small, the VDOS of the 
SWNT deviates from that of the graphite due to the pe- 
riodical boundary along the circumferential direction and 
the curvature effect, and the amplitude of the deviation 
is almost inversely proportion to n. As a result, when (n, 
0) and (2n, 0) tubes connect with each other, the differ- 
ence between the VDOSs of the two tubes is obvious if 
n is small, and vice versa, the difference is not so obvi- 
ous if n is large. Recalling the fact that the rectification 
can be related with the mismatching between the VDOSs 
around the connecting region, we can conclude that the 
rectification of the IMJ is not obvious when n is large. 

Next, we investigate the length dependence of heat 
conduction. From Fig. 6a, it can be seen that the total 
heat flux decreases as the system length increases from 
Lq to 2Lq. However, the thermal conductivity, k — 2 'sAT > 
increases with the increasing of the system length. This 
result agrees with the previous studies. Moreover, there 
still exists asymmetry in the heat flux, which is clearly 
illustrated in Fig. 6b. Where the rectification at different 
| AT | is plotted for different system lengths. Obviously, 
the increase of system length weakens rectification. It 
can be understood as follows. As the system length in- 
creases, long-wave phonons modes contribute more to the 
heat conduction. This kind of heat conduction is sym- 
metric because long-wave phonon modes are hardly scat- 



tered by the junction or any local defects. Moreover, the 
asymmetric thermal conductance is mainly controlled by 
the asymmetric interface, i.e., the connection part, which 
does not change when system length increases. As a re- 
sult, although the absolute value does not change, the 
rectification appears to be smaller. 

Based on above study, we can conclude that changes in 
radial and axial direction have different influence on rec- 
tification. The increase of structural asymmetry around 
the interface can increase the rectification. In the follow- 
ing section we will discuss other alternatives to improve 
the rectification. 



IV. IMPROVEMENT 

First, it is possible to increase the difference in force 
constants around the interface by applying external 
stress. Because the Young's modulus has been proved to 
be almost the same in different tubes using the Tersoff- 
Brenner potential^ it can be expected that the (n, 0) 
tube and (2ro, 0) tube will give different response to the 
same external stress. Thus it is possible that one can 
use this property to fabricate devices. Here we consider 
two kinds of stresses, tensile and torsional stresses. It 
has been proven in the framework of tight-binding*- that 
these two stresses can change the electronic structure of 
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SWCNTs greatly. Here the (8, 0)/(16, 0) SWCNT IMJ 
(to = 4) is taken as an example. The original length is 

The stresses considered in this work are all small 
enough to ensure the system does not undergo any struc- 
tural transformations^ 4 - or defects^ For tensile stress, 
we elongate the IMJ along the axial direction by 3%. 
For torsional stress, one head of IMJ is rotated relative 
to the other by 30 ° around the axis. The results are 
shown in Fig. 7. 
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FIG. 7: (Color online). The thermal transports of (8, 0)/(16, 
0) SWCNT IMJs under different stresses, (a) The heat flux 
for different stresses, (b) The thermal rectification for differ- 
ent stresses, m = 4. Error bars are also plotted. 

It can be seen that the torsional deformation can 
change neither the absolute value of heat flux nor the 
asymmetric behavior largely. However, the tensile stress 
can greatly change the heat transport. In the elongated 
structure, the heat flux becomes smaller than that of un- 
deformed structure, but the rectification becomes larger. 

Above results are reasonable. First of all, the defor- 
mation of the thicker tube is smaller than that of thin- 
ner ones under same external stress. This means that 
when thick and thin tubes are connected with each other 
to form IMJ, they will appear different axial deforma- 
tion under the same stress, which will further increase 
the structural difference between the two tubes. Sec- 
ondly, the heat transport in carbon nanotubes are mainly 
contributed by the four acoustic phonon modes. These 
phonon modes have weak dependence on the local struc- 
ture, but they can be affected by the shape of the cylinder 
surface. When tensile stress is applied to the SWCNTs, 



the SWCNTs will shrink in the radial direction and elon- 
gate in the axial direction. This means their cylinder 
surface will change, so that the acoustic phonon modes 
will also be changed. But when SWCNTs are twisted, 
the cylinder surface almost does not change, so that the 
thermal conductance almost does not change. 

Comparing Fig. 7a and 3a, it can be found that the 
heat flux decreases in elongated structure. The reason 
is straightforward. As mentioned above, the thermal 
conductivity can be expressed as k = ^ l\ q C q (u>\) v\ q , 

i 

For the low 
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where C q (uj) - k B ^ j-? ) [cxp{h u/k B T)-\F 
frequency (hw <C ksT) modes, which are most important 
for the thermal transport, C q (to) ~ ks- Furthermore, l\ q 
and v\ q are the increasing functions of the force constants 
k, thus the thermal conductivity k is also an increasing 
function of the force constants k. In fact, some numerical 
simulations show that k oc k 2 in the weak coupled lat- 
tice models^i 4 ^ When the system is elongated, the bonds 
along the axis are also elongated. Then according to 
the Abell-Tersoff formalism of the Tersoff-Brenner po- 
tential, the corresponding force constants will decrease. 
Finally, the thermal conductivity decreases in the elon- 
gated structure. 

To elongate an IMJ is not very difficult for the state- 
of-the-art experiments, so we believe that tensile stress 
can be a possible method to improve the rectification. 

Next, we will try to modify the structure to change 
rectification. 

It is suggested^ that a suitable on-site potential can 
induce large rectification in two-segment model. So it is 
quite interesting that what will happen if some period- 
ical potential is introduced into the SWCNT IMJs. On 
the other hand, a nanotube with sufficiently large diam- 
eter can be filled with spherical Cqq fullerene molecules 
to build up a new hybrid structure referred to as a 
'peapod' , 48 ' 49 with spherical fullerenes representing peas 
and the carbon nanotube representing a pod. The ther- 
mal conductivity of the peapod structure has been inves- 
tigated by using MD recently— So it is possible to com- 
bine these two structures, i.e., IMJs and peapod, into a 
new structure. 
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FIG. 8: (Color online). The (8, 0)/C 6 o@(16, 0) structure. 
The number of pent agon- heptagon defects is m = 4. The 
regions marked as T are fixed in MD process. The regions 
marked as 'II' are put in the heat baths. 

The structure (see Fig. 8) contains two parts, i.e., 
a segment of (n, 0) SWCNT and a segment of (2n, 0) 
SWCNT with some Ceo balls inserted into its center, 
which is typically called as SWCNT peapod structure. 
The lengths of these two segments are almost the same. 
They are connected by to pairs of pentagon-heptagon de- 
fects. Because the closest distances between the adjacent 
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Ceo balls are almost same, 10.05 A, we can consider that 
they exert an external periodic potential on the outside 
SWCNT. Here, the interaction between two atoms on Ceo 
balls and outside SWCNT is modeled by Lennard-Jones 
potential: 



V (r) = 4e 



<J\ 6 

r / 



r 



(3) 



with parameters, e = 2.964 meV and a = 3.407 A. By 
default, the total length of the structure is Lq, which 
permits 10 Ceo balls in the (2n, 0) segment. 

In this complex structure, the two heads are closed to 
ensure the Cgo balls not running out in the MD calcu- 
lations. And several periods of two heads (presented by 
blue color and marked as region T in Fig. 8) are fixed 
in the MD simulation. Then two periods of each ends 
(emphasized by orange color and marked as region 'II' in 
Fig. 8) are put in the Nose-Hoover thermostat^ 

In Fig. 9, the heat transports of (8, 0)/C 60 @(16, 0) 
structure with different lengths are shown. 




0.4 0.6 

I AT|/(r) 

FIG. 9: (Color online). The thermal transports of (8, 
0)/C6o@(16, 0) structures with different length, (a) The heat 
flux versus temperature difference, (b) The thermal rectifica- 
tion versus temperature difference, m — 4. Lo is defined in 
text. Error bars are also plotted. 

Comparing Fig. 9a with Fig. 6a, one can find that 
the heat flux of peapod structure is smaller than that of 
the empty IM J, which seems to be in contradict with the 
result in Ref. l50l . But it should be pointed out that in our 
structure the peapod region is closed at the ends. In fact, 



this prevents the possible mass transport (via fullerenes) 
which in turn reduces the assistance transport effect of 
Ceo balls. Furthermore, because of the coupling between 
the phonon modes (mainly the radial breathing modes) 
of carbon nanotube (CNT) and Cgo balls, the conduction 
phonon will be scattered and the thermal conductivity of 
outside CNT will be reduced. 

On the other hand, we can notice that the length of 
error bars in Fig. 9b is much larger than that in Fig. 
6b. This means that the existence of Cgo balls increases 
the fluctuations in the heat flux. In Ref. 50j, it is found 
that the temperature fluctuations of peapod structure 
are larger than that of empty nanotube. It is explained 
by the weak van der Waals forces between fullerenes and 
the outside CNT. 

Another result is that the existence of Ceo balls can 
hardly change the asymmetric behavior of heat flux. This 
is somehow disappointing for us, because it is natural 
to expect that the periodical potential introduced by 
the Ceo balls can enhance the thermal rectification ac- 
cording to the results based on lattice models. How- 
ever, our numerical results are also reasonable because of 
following facts. The periodical potential introduced by 
the Ceo balls is intrinsically van der Waals potential, so 
its strength is much weaker than the binding potential 
among the carbon atoms in SWCNT. The period of the 
periodical potential introduced by Cgo balls is about 10 
A, which is too large compared with the C-C bond length 
in SWCNT, thus only low frequency vibrational phonons 
are affected. As has been mentioned, the asymmetric 
heat flux is controlled by the optical phonon modes, so 
the change of low frequency vibrational modes can hardly 
change the asymmetric behavior. 

According to our result, if one wants to build a thermal 
rectifier by applying periodical external potential, van 
der Waals force is not a good choice. 



V. CONCLUSIONS 

In this work, heat rectification in (n, 0)/(2n, 0) IMJs 
has been studied. It is found that the rectification de- 
pends weakly on the detailed structure of connection 
part, but depends strongly on the temperature gradient. 
The rectification increases as the temperature gradient 
increases. 

We have also studied the dependencies of the rectifi- 
cation on tube radius and IMJ length. We found that 
the maximum rectification appears in short (8, 0)/(16, 
0) IMJ. 

Moreover, we have found that the tensile stress can 
increase the rectification, while torsional stress can not. 
It is due to the different effects of the stresses on IMJ 
structure. It can be believed that the tensile stress will 
play an important role in building of devices with highly 
asymmetric transport behavior. 

The 'peapod' structure is combined with IMJ structure 
to investigate the effect of periodical potential on the 
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rectification. Our numerical result shows that the new Agreement No. POD0410553. 
structure can hardly change the rectification. 

The study may shed lights in further experimental in- 
vestigation in this field. 
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